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ABSTRACT 



We present new results from a significant extension of our previous high angular 
resolution ((X'3 ~ 40 AU) Submillimeter Array survey of the 340 GHz (880 /xm) thermal 
continuum emission from dusty circumstellar disks in the ~lMyr-old Ophiuchus star- 
forming region. An expanded sample is constructed to probe disk structures that emit 
significantly lower millimeter luminosities (hence dust masses), down to the median 
£/^ \ value for T Tauri stars. Using a Monte Carlo radiative transfer code, the millimeter 

visibilities and broadband spectral energy distribution for each disk are simultaneously 
reproduced with a two-dimensional parametric model for a viscous accretion disk that 
has a surface density profile S oc {R/Rc)^ 1 exp [— (R/Rc) 2 ^ 1 ]. We find wide ranges of 
characteristic radii (R c = 14-198 AU) and disk masses (M d = 0.004-0.143 M ), but a 
narrow distribution of surface density gradients (7 = 0.4-1.1) that is consistent with a 
uniform value (7) = 0.9 ± 0.2 and independent of mass (or millimeter luminosity). In 
this sample, we find a correlation between the disk luminosity/mass and characteristic 
radius, such that fainter disks are both smaller and less massive. We suggest that 

o 



this relationship is an imprint of the initial conditions inherited by the disks at their 
formation epoch, compare their angular momenta with those of molecular cloud cores, 
and speculate on how future observations can help constrain the distribution of viscous 
evolution timescales. No other correlations between disk and star properties are found. 
The inferred disk structures are briefly compared with theoretical models for giant planet 
formation, although resolution limitations do not permit us to directly comment on 
material inside R ~ 20 AU. However, there is some compelling evidence for the evolution 
of dust in the planet formation region: 4/17 disks in the sample show resolved regions of 
significantly reduced millimeter optical depths within ~20-40 AU of their central stars. 

Subject headings: accretion, accretion disks — circumstellar matter — planetary sys- 
tems: protoplanetary disks — solar system: formation — stars: pre-main-sequence 
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Introduction 



Direct observations of the reservoirs of planet-building material — the disks around young stars 
— should play a critical role in developing theoretical models of planet formation. Ultimately, the 
goal of those models is to elucidate the physical processes involved in making planets by incorpo- 
rating the observed properties of circumstellar disks and successfully reproducing the demographic 
characteristics of the p lanets in our Solar System and those around other stars (e.g-. llda &: Linll2004l : 



Mordasini et al. 



200 



Regardless of the favored mechanism, two basic requirements for forming 
planets in a circumstellar disk must be satisfied: there must be enough material (gas and dust) in 
the right locations and a suffic ient amount of time for the formation mechanism to operate (e.g., 
Pollack et al.lll996l ; iBossI 119971 ). The former criterion amounts to a density threshold, suggesting 
that observational constraints on the spatial distribution of mass in young circumstellar disks are 
fundamental in constructing an empirical foundation for planet formation models. Moreover, those 
same estimates of disk densities can be used to help characterize th e viscous accretion process that 
determines how disk structures evolve (e.g.. iHartmann et al.lll998l ). 



There are some significant observational obstacles to direct measurements of disk densities. 
Most of the material in these disks is "dark," composed of cold molecular hydrogen and not readily 
detectable. Their mass contents must be inferred from trace species, particularly from the dust 
grains that dominate the disk opacity. At radio wavelengths, the thermal continuum emission 
from these dus t grains is optically thin and therefore provides an unique probe of mass in the 
disk midplane (JBeckwith et al.lll990i ). If that emission can be resolved, properly interpreted with 
radiative transfer calculations, and assigned some nominal gas-to-dust mass ratio, it can be used to 
reconstruct the spatial distribution of disk densities. Interferometric observations o f dust and trace 



gas species have helped quantify d e nsities in the outer parts of circumstel lar disks (JKitamura et al. 



2002 



Andrews &: Williams 



2007a 



Hughes et al. 



2008 



presence of vertical temperature and density gradients (jDartois et al.1 



Isella et al. 



2009), clearly establishe d the 



2003 



Pietu ct al. 



even identified disks with large central cavities that exhibit very little du st emission (jPietu et al. 



2006 



Hughes et al. 



2007 



2009 



Brown et al.ll200Sl . 120091 : llsella et alJl2010h . 



2007) . and 



In a previous study (j Andrews et al.ll2009l . hereafter Paper I), we presented the initial part of a 
high angular resolution (0 / /3 ~ 40 AU) survey of the 880 /an thermal dust emission from protoplan- 
etary disks in the nearby Ophiuchus star-forming region. To take advantage of these observations 
that are sensitive to the mass content inside the planet formation zone (R < 40 AU), we devel- 
oped a radiative transfer modeling toolkit to extract the two-dimensional temperature and density 
structures of these disks through a combined fit of their millimeter continuum visibilities and broad- 
band spectral energy distributions (SEDs). Through that modeling effort, we concluded that the 
disks had densities comparable to those expected for the outer parts of the primordial disk around 
the Sun and very similar radial density gradients. The inferred density gradients indicate that 
whatever mechanism is responsible for generating viscosity in these disks acts with a nearly linear 
radial distribution (i.e., the viscosities vary as v ex R? , with a median 7 ~ 0.9). Moreover, we 
resolved large regions (R ~ 20-40 AU) with very little millimeter emission in the centers of several 



disks, and speculated that these may be the signposts of young planetary systems. However, that 
initial sample was biased towards the targets that are exceptionally bright at millimeter wave- 
leng ths, and therefore considera bly more massive than the typical disk in the Ophiuchus region 
(see I Andrews fc Williams! l2007bl ) . In this article, we extend the survey to double the sample and 
probe disk structures with millimeter luminosities down to the median value for ~lMyr-old T 
Tauri stars. In §2 we describe the sample selection criteria, new observations, and data calibration. 
A brief review of the radiative transfer calculations is provided in §3, with the modeling results for 
this extension of the sample highlighted in §4. In §5, we synthesize the disk structure constraints 
for the full sample in the contexts of the viscous evolution process and the prospects for planet 
formation. The key conclusions from the survey are summarized in §6. 



2. Sample Selection, Observations, and Data Reduction 

The sample of disk targets for this survey was selected primarily for observational and analyt- 
ical convenience. To ensure that the SMA data would be sufficiently sensitive to probe emission 
over a large range of spatial scales, targets were required to have integrated 880 /im flux densities 
larger than ^75mJy. Candi dates were selected based on the single-dish photometry compiled by 



Andrews &: Williams! (|2007bl ). When no 880 um data were available, we relied on the 1. 3 mm flux 
densities provided by previous surveys ( Andre fe Montmerldll994 : lNiirnberger et al. 19981) and con- 
servat ively scaled up by the square of the wavelength ratio (~2.2; see I Andrews fe Williams! l2005bl . 
2007b! ) . After this initial cut, we excluded targets that lacked sufficient information about their 
central stars. Without stellar temperatures and luminosities, it is not possbile to interpret the ob- 
served millimeter data in detail (see §3). In practice, this second criterion amounted to an extinction 
threshold Ay < 15, as spectral c lassifications and luminosity estimates are rare and uncertain fo r 



more deeply embedded sources (JLuhman &; Rieke 



1999 



Wilking et al. 



2005 



Furlanet al J 12003 ). 



The sample selection-space set by these criteria is shown in Figure (H with red points marking the 
selected target disks. The 880 /im flux densities for those points correspond to the values derived 
from the SMA observations described here or in Paper I. One selected target, WL 18, falls below the 
designated millimeter emission threshold: previous 1.3 mm data led us to expect a larger 880 /im 
flux density than was observed. One other target, RX J1633.9— 2422, meets t he selection crite ria 
but was not observed: its millimeter flux density was only recently published (jCieza et al.ll2010i ). 



These combined criteria yielded 17 disk targets for the sample. The 880 /im flux density 
criterion corresponds to the median yalue fo und in both the Ophiuchus and Taurus star-forming 
regions ([Andrews fc Williams! l2005bl . l2007bl ). Therefore, the sample defined here fully spans the 
upper half of the millimeter continuum luminosity distribution — or equivalently the circumstellar 
dust mass distribution — for the disks around ~1 Myr-old stars. For the typical assumptions 
in converting this emission to total disk masses, this limit corresponds to a few Jupiter masses 
(~0.004M Q ) of gas and dust, suggesting that the sample is representative of the disk population 
that might eventually be able to produce giant planets. Moreover, the sample includes targets that 



essentially span the full range of T Tauri star properties, including spectral types from M4 to G3, 
stellar masses of 0.3-2.0 M , and accretion rates of ~10 _9 -10 _7 Mq yr _1 . 

To complete this sample and expand on the initial 9 disks discussed in Paper I, 8 additional tar- 
gets were observed with the very extended (V; 8-509 m baselines) and compact ( C: 6-70 m baselines) 



configurations of the Submillimeter Array interferometer (SMA; iHo et al.ll2004l ) in 2009 March and 
May. A journal of these SMA observations is provided in Table [TJ The SMA double sideband 
receivers were tuned to a local oscillator (LO) frequency of 340.755 GHz (880 /urn). Each sideband 
was divided into 24 partially overlapping 104 MHz chunks centered ±5 GHz from the LO frequency. 
The central chunk in the upper sideband was sampled at a factor of 4 higher spectral resolution than 
the others, in an effort to observe the CO J=3— 2 transition (345.796 GHz) in 0.70 km s _1 channels. 
The observing sequence interleaved disk targets with nearby quasars, J1625— 254 and J1626— 298, 
in an alternating pattern with a total cycle time of ~10-15 minutes. When the targets were at 
low elevations (<20°), planets (Uranus, Saturn), satellites (Titan, Callisto), and bright quasars 
(3C 454.3, 3C 279) were observed as bandpass and absolute flux calibrators depending on their 
availability and the array configuration. The observing conditions were generally very good, with 
atmospheric opacities <0.1 at 225 GHz (corresponding to <2.0mm of precipitable water vapor). 

The data were edited and calibrated as in Paper I using the MIR software package. The bandpass 
response was calibrated with observations of a bright planet or quasar, and broadband continuum 
channels in each sideband were generated by averaging the central 82 MHz in all of the chunks 
except the one reserved for the CO J=3— 2 line. The visibility amplitude scale was set based 
on observations of planets or satellites and routinely-monitored quasars: the typical systematic 
uncertainty in the absolute flux scale is ~10%. The antenna-based complex gain response of the 
system as a function of time was determined with reference to J1625— 254, which lies only ~1° from 
the target disks. The other quasar in the observing cycle provides a check on the quality of the 
phase transfer in the gain calibration process. The millimeter "seeing" generated by atmospheric 
phase noise and any small baseline errors is small, 0.1-0 / /2. After combining all of the data for 
each target, the standard tasks of Fourier inverting the visibilities, deconvolution with the CLEAN 
algorithm, and restoration with a synthesized beam were conducted with the MIRIAD software 
package. High resolution maps of the continuum emission were created with a Briggs robust = 
0.2-0.7 weighting scheme for the visibilities, and maps of the CO J=3— 2 line emission were made 
with natural weighting for the compact array data only. The relevant data properties from these 
synthesized maps are compiled in Table [2j 

The synthesized continuum maps for these targets are featured in Figure [2j Each high angular 
resolution map covers A" on a side, corresponding to 500 AU at the adopted distance of 125 pc to 
the O phiuchus clouds (Ide Geus et al.lll989i : iKnude fc HogHl998l : lLombardi et al.ll2008l : lLoinard et al 



2008). A centroid position and initial estimate of the viewing geometry of the disk — characterized 
by the inclination (z) and major axis position angle (PA) — were determined by fitting the visibilities 
with an elliptical Gaussian brightness distribution. Because the continuum emission from the SR 
24 system is not centrally peaked, its centroid position and viewing geometry were estimated by 



inspection of the lower resolution (compact SMA configuration) data only. In all cases, the CO line 
emission from the disks is significantly contaminated by the local molecular cloud environment, 
and will not be discussed further. 



3. Modeling the Disk Structures 

To interpret these observations in the context of the disk structures — physical conditions 
and key size scales — present in this sample, we followed the modeling formalism introduced and 
discussed in detail in Paper I. To briefly summarize that procedure, we first adopt a parametric 
prescription for a two-dimension al, flared density structure based on a simple model f or the viscous 
evolution of an accretion disk (JLvnden-Bell &; Pringld 11974 ; lHartmann et al.lll998l ). The model 
assumes an anomalous disk viscosity that varies with radius as v oc W and has a surface density 
profile 



(2-7) 
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exp 
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R r 



(1) 



where R c is a characteristic scaling radius and Md is the disk mass. At a given radius, that 
column density is vertically distributed so that it falls off with altitude above the midplane like a 
Gaussian with scale- height H oc R 1+ ^. Once a density structure has been specified with a set of 
five parameters, {Md, 7, R c , H c , ip} (where H c is the scale-height at the characteristic radius), it is 
populated with a spatially homogeneous size distributio n of dust grains. As described in Paper I, we 
consi der dust with an interstellar medium composition (JDraine fc Ledll984l : IWeingartner fe Draine 
200 ll ) and a power-law distribution of sizes (with an index of -3.5) between 0.005 /xm and 1 mm. 
This dust structure is then irradiated by a central star of fixed temperature and luminosity to 
determine an internally-consistent disk temperature structure, u sing the two-dimensional Monte 
Carlo radiative transfer code RADMC (JDullemond &: Dominikll2004l ). The results of the Monte Carlo 
simulation are then coupled with a post-processing raytracing code to compute a synthetic dataset 
for a given viewing geometry, consisting of both a broadband spectral energy distribution (SED) 
and a set of millimeter continuum visibilities that sample the Fourier plane in the same way as 
the SMA data. These synthetic data are simultaneously compared with the observations, and the 
process is iterated until the fit converges on a minimum joint x 2 value (see Paper I for details). 

The stellar properties used to generate the Kurucz spectra that irradiate the disk structures 
are listed in Table [3l based on both literature measurements and matches to the optical/near- 
infrared SED. Because the spatio-kinematic information from the CO line emission in these disks 
is sufficiently contaminated by the local cloud material, estimates of the disk inclinations and 
position angles were determined solely from the elliptical Gaussian fits to the continuum emission 
described in §2. Unfortunately, the relatively faint continuum signals from some of these targets 
leads to considerable uncertainty in these viewing geometry estimates. These are still the best values 
available, but we caution that some systematic uncertainty will remain for the density parameters 
until better emission line data are available for these targets. With such a limited amount of 



spatially resolved information about the two components in the DoAr 24 E system, a detailed 
modeling of this source is beyond the scope of this study. Some important modifications to the 
modeling process were made for the SR 24 S disk, with the details provided in §4.2. 



4. Results 

4.1. Disk Structures 

The model parameter values that best reproduce the observations for the full sample are 
compiled in Tables 0] and [5] for the disks with continuous emission distributions and those with 
central emission cavities, respectively. For the continuous disks, we include the (fixed) inner disk 
radii inferred from a simple sublimation argument (see Paper I; Tabled column 7). Columns 7 and 
8 in Table [5] list the estimated sizes (-R C av) and density contrasts (8 cav ) used to account for the disks 
with central emission cavities. Both tables also include the adopted values for the disk inclination 
and position angle, as well as the reduced \ 2 statistics for the fits to the SED and visibility datasets 
separately. The distributions of these values for each parameter are shown together in Figure [3j 
hatched regions mark contributions from the disks with central cavities (around SR 24 S, SR 21, 
DoAr 44, and WSB 60). In Figure 0J the new observations presented here are directly compared 
with the synthetic datasets generated from these best-fit models (see Paper I for the other sample 
disks). From left to right, we display the observed SMA millimeter continuum image (as in Figure 
[2]), the synthesized model image, the imaged residuals, the broadband SED, and the elliptically 
averaged millimeter visibility profile (see Paper I for details). The latter two panels have the best- 
fit model behavior overlaid in red, and the SED panel also shows the input stellar spectrum as a 
dashed blue curve. Because of its large central emission cavity, the modeling results for the SR 24 
S disk are shown separately in Figure with inset images synthesized at higher angular resolution 
to facilitate a more detailed comparison. 

The modeling uncertainties and relationships between the observations and model parameters 
were already discussed in detail in Paper I, so we will not repeat that effort here. Instead, we 
focus on how this expanded sample enables new constraints on the spatial distribution of mass in 
~1 Myr-old circumstellar disks. The inferred surface density profiles for the full sample are shown 
together in Figure [6l with the targets presented for the first time here highlighted in color. The 
disks with large central emission cavities are shown in a separate panel. The light gray band inside 
20 AU marks the survey resolution limit, while the dark gray boxes are reference points representing 
the minimal surface densities (with uncertainties) and presumed feeding zones for Satu rn, Uranus, 



and Neptune in the canonical model for the primordial solar disk (j Weidenschillingj 1 1 9771 ) . Although 
the disks in this sample exhibit a wide range of masses (M^; Figure EK) and characteristic radii 
(R c ; Figure [3b), we find a very narrow distribution of values for 7, the radial gradient of the surface 
density profile. As shown in Figure 03, that distribution is roughly normal with a peak at 7 = 0.9 
and a standard deviation of 0.2. Since that width is comparable to the modeling uncertainties on 



7 for an individual source (~0.2-0.3; Paper I), the modeling results are consistent with a uniform 7 
value being representative of the entire sample^] Moreover, by extending the survey to cover much 
fainter targets, we have confirmed that the shape of the surface density profile does not significantly 
change over a wide range of millimeter luminosities (or disk masses) . 

Despite the similar radial density gradients in these disks, the observations clearly show a wide 
variety of millimeter emission morphologies. Some of that diversity is related to how the disks 
are heated by their central stars, a natural outcome of the assortment of stellar properties and 
vertical distributions of dust present in the sample (see Paper I for details). However, the density 
profiles plotted in Figure [6] demonstrate that a significant range of characteristic radii and masses 
are also partly responsible. A close examination of those S profiles reveals that the brighter disks 
at millimeter wavelengths tend to have both higher masses and larger characteristic radii. The 
former is no surprise, since the vast majority of the disk volume at these wavelengths is optically 
thin, but there is no a priori reason to expect the disk sizes to be correspondingly smaller for fainter 
sources. This relationship is not a modeling artifact, as there is an analogous empirical correlation 
between the brightness of the millimeter emission and how well that emission is resolved. Figure [7] 
demonstrates this explicitly by comparing the average 880 /mi visibility profiles for disks that are 
brighter (red) or fainter (blue) than 0.5 Jy. These profiles were calculated by averaging the visibilities 
for each subset of disks into annular bins, after deprojecting the visibilities for individual disks 
according to their viewing geometries and normalizing their real (correlated) fluxes by dividing 
off their integrated flux densities. This comparison confirms that the millimeter emission from 
the brighter disks in this sample is more resolved by our SMA data, exhibiting substantially less 
correlated flux on essentially all spatial scales. Although the physical origins of this empirical 
relationship are not obvious, some speculations about its significance for understanding the viscous 
evolution process in these disks are presented in §5. 

In an effort to identify any trends among the disk structures or the characteristics of their 
stellar hosts, we performed a principle component analysis on a subset of such properties for the 
12 disks with continuous density distributions in this sample. The analysis included the 5 free 
disk structure parameters, {M^, 7, R c , -ffiotb V 7 }; the stellar properties {T e ff , L*}, and the accretion 
rates, M* (see TablesEHH and[6|). Due to the additional uncertainties of the stellar evolution models 
used to infer them, we did not include stellar ages, masses, or radii directly (although those values 
are well-correlated with L* and T e g). As alluded to above, we identified one statistically significant 
correlation (3.3 a; correlation coefficient of 0.85) between M^ and R c , which will be discussed in 
detail in §5.1. The first eigenvector from the principle component analysis, accounting for 40% 
of the variance in the data, is dominated by a positive trend relating R c , Mj, L*, and M*. The 
trend indicates that the more massive disks in this sample are larger and orbiting more luminous - 
perhaps younger - stars that are accreting disk material at higher rates. While these relationships 



The only major outlier, the disk around AS 209, has a particularly uncertain 7 estimate due to a poorly-sampled 
infrared SED (see Paper I for details). 



may hint at crucial information related to the viscous evolution process, a larger sample will be 
required to make any definitive conclusions. 



4.2. Commentary on Individual Disks 

Elias 24 — This heavily-reddened classical T Tauri star is one of the brightest millimeter con- 
tinuum sources in Ophiuchus. In fact, it was initially excluded from the sample out of concern 
that the bright emission signaled contamination from an extended envelope. However, there 
is no evidence for the spatial filtering of such large-scale emission; the SMA flux densities at 
both 880 //m and 1.3 mm are i n excellent agreement with lower resolution single-dish photometry 



( Andrews fc Williams 



visible (Wil king et al. 



lope (JBontemps et al. 



2007al Jbl; lAndre Sz Montmerld 11994 ). Moreover, the central star is optically 



2001 



2005) and the infrared SEP shape is incomp atible with a substantial enve- 



I Evans et al. 



2003 



Barsonv et al.ll2005l ). The 880 /mi continuum data 



presented here features a bright central core and a fainter, extended, and apparently asymmetric 
emission halo on larger scales. Although the best-fit model for this source is able to reproduce the 
SED and visibilities rather well, any axisymmetric model will underpredict the extended emission 
to the east of the disk center at the ~3a level (see Figure H|) . The origins of that extension remain 
unclear. 



SR 24 S — The hierarchical triple s ystem SR 24 is composed of this K2 s tar and a close binary 



pair (Of! 2 « 25 AU; ISimon et al.lll995l ) located ~5" (625 AU) to the north (IReipurth & Zinnecker 



1993). Both SR 24 S and the SR 24 N binary exhibit excess emission from warm dust disks, 
bright Ha lines indicative of substantial accretion flows, and extended emission from CO low- 
energy rotational tra nsitions. However, all of the continuum emission at millimeter way elengths is 
produced by SR 24 S d Andrews fc Williamsll2005al : IPatience et al.ll2008l : llsella et al1l2009h . The high 
angular resolution inset image of the SR 24 S disk in Figure [2] reveals a resolved central emission 
cavity with an apparent brightness enhancement to the northeast. The origins of this ring-like 
emission morphology are unclear, but could perhaps be generated by abrupt emissivity variations 
due to particle growth, a dramatic dissipation process driven by high-ener gy radiation from th e 
star, or even tidal interactions with companion objects (in this last case, see lMayama et al.ll2010l ). 



Regardless of the cause, the model described in §3 was adjusted to account for the observed 
disk morphology. As detailed in Paper I for similar cases, we adopted the simple modification 
of artificially decreasing the surface densities inside a radius R ca _ v by a factor (?„ a nd then at- 



tempte d to fit the SMA visibilities and the component-resolved optical ( Wilkin g et al.1 120051 ) , in 



frared (I.Tensen fc Mathiei]lll997l : ICutri et al.ll2003l : Evans et al.ll2003l : lMcCabe et al. 



3), 



2006 , and radio 



photometry. However, much like the case of the DoAr 44 disk described in Paper I, it is difficult 
to simultaneously account for the significant infrared excess and the lac k of millimeter c ontin uum 
emission near the SR 24 S stellar position. Following the approach of lEspaillat et al.l (|2007l ) for 
other "pre-transitional" disks, we proceeded by first finding a good match to the millimeter visi- 



bilities, and then artificially increasing the surface densities near the star until sufficient infrared 
emission was produced to match the SED (without affecting the millimeter emission). A satisfac- 
tory result was achieved with a surface density profile scaled down by 5 CSV ss 0.05 inside R = 2 AU, 
a factor of 20 lower than for a continuous disk, but ~500x higher than for the region between 2 AU 
and -R cav ~ 32 AU. Obviously these adjustments are artificial and non-unique; a modeling effort 
focused on a more robust exploration of the detailed inner disk structures will be treated elsewhere. 



SR 4 — This star harbors a disk with a compact, centrally-peaked millimeter emission morphology. 
The low foreground extinction allows for a particularly good characterization of the stellar lumi 



nosity from optical photometry, despite the low-level variability at those wavelengths (IHerbst et al 



1994 ). The infrared SED s hows a relatively small excess at sho rt wavelengths before a pronounced 



flattening from ~8-24 / um (jCutri et al. 



2003 



Evans et al.ll2003l ). Unfortunately there is little addi- 



tional information in the far-infrared to better constrain the SEP m orphology, as the Spitzer 70 /xm 
images are contaminated by local nebulosity ([Padgett et al.ll2008l ). This lack of data between the 
mid-infrared and millimeter regions of the SED results in only a relatively crude constraint on the 
vertical structure parameters (and thus temperatures) in this case. 



SR 13 — Another hierarchi cal triple system ( see ISchaefer et al.ll2006l ). SR 13 includes a close 
binary p air (13 mas ss 1.6 AU; ISimon et al.lll995l ) and a tertiary companion located ~0 / .'4 (50 AU) 
due east (JGhez et al.lll993l ). Lacking sufficient component-resolved photometry for this system, we 
made a preliminary model ignoring the tertiary and assuming the primary is a s ingle star. The 
SED was constructed from the composite op tical measurements of IHerbst et al. (119941 ) . infrared 
data from 2MASS and the Spitzer c2d survey (ICutri et al.ll2003l: lEvarts et al.l l2003). and single-dish 



radio observations (jAndre &: Montmerld [1994J ; lAndrews &: Williams! l2007bl ) . Despite this simplifi- 
cation, the SMA data can potentially resolve another disk around the tertiary. There is a modest 
eastern extension in the emission map, and a potential null near 400 kA in the deprojected visibility 
profile (Figure |4]). These features could be produced either by a cavity in a circum-system disk 
or by a marginally resolved additional disk around the spatially offset tertiary. With the available 
sensitivity, it is difficult to differentiate these scenarios. Needless to say, the model disk structure 
inferred in this case should be treated with appropriate caution. 



WSB 52 — This classical T Tauri star with spectral type Ml exhibits a compact millimeter 
continuum emission distribution typical of the fainter disks in t he sample. The SED used in 



the model fits was construct ed from the red-o ptical photom etry of IWilking et al.l (120051). infrared 



measu rements from 2MASS (ICutri et al.ll2003 ). Svitzer c 2d ([Evans et al .1120031) . and Padgett et al. 
( 20081 ). and a 1.3 mm single-dish flux measurement from lStanke et al.1 ( 2006 ). 



DoAr 33 — The infrared excess emission around this K4 star is remarkably f aint compared 



to the typical young disk. We constructed a SED from red-optical photometry (IWilking et al. 



10 



20051 ), the 2MASS an d Spitzer c2d programs (ICutri et al.l 12003: Evans et aljj200j ), and single-dish 
radio measurements ([Andre fc Montmerld Il994l ; lAndrews &: Williamsl l2007bl ) . Only a very weak 
continuum excess is present shortward of ~8 jttm. At longer wavelengths where the dust excess is 
brighter, the spectrum is very steep and blue. Reproducing this infrared SED shape with our models 
required the u se of a very f l at ver tical distribution of dust, resulting in comparatively cold midplane 
temperatures. ICieza et al.1 (120101 ) argue that these infrared colors are indicative of significant dust 
evolution. The flat structure inferred here could be interpreted in that context as evidence for 
advanced dust sedimenation to the midplane, or perhaps the sign of extensive shadowing of the 
outer disk due to some perturbed inner disk structure. In any case, the coupling of the rare infrared 
SED, bright millimeter emission, and weak signatures of accretion (jCieza et al.ll2010i ) suggest that 
this disk is worthy of further scrutiny. 



WL 18 — This star harbors the faintest disk in the sample, due in part to its small size (R c = 
14 AU). Its visibility profile shows that the millimeter emission is only marginally resolved on the 
longest SMA baselines. The SED compiled here is spars e due to high for e ground extinction , with 



infrared data from 2MASS and the Spitzer c2 d project (ICutri et al. 



a single-dish 1.3 mm flux density measured bv lMotte et al.l (J1998I ). That same extinction makes an 



2003 



Evans et al.ll2003l ). and 



estimate of the underlying stellar luminosity difficult, which in turn contributes an uncertainty to 
the thermal structure of the disk. The apparently low luminosity of the star relative to others with 
the same spectral type (K7) leads to an uncomfortably large age estimate (11 Myr). 



DoAr 24 E — This spectral t ype G6 weak-lined T T au ri star has a companion roughly 2" to 
the southeast (JGhez et al.lll993l : iReipurth fc Zinneckerlll993r) that becomes sub stantially brighter 



than the primary at wavelengths longer than ~3/^m (e.g., iMcCabe et al.ll2006l ). The unresolved 



composite SED has been attributed to the photosphere of the primary and dust around this infrared 



Barsony et al. 



2005 



companion. However, reso l ved measurements a t near- and mid-infrared wavelengths (jPrato et al 

200; 



McCabe et al.1 120061 ) indicate that the primary does show some very 



weak infrared excess at least at ~10//m. Regardless, most of the unresolved millimeter emission, 
and therefore the c i rcumstellar dust mass , was as sumed to be generated by the infrared companion 
(jMotte et al.l Il998l : lAndrews fc Williamsl l2007bl ) . The resolved 880 /jm image of this system in 
Figure [2] clearly demonstrates that this is not the case: each stellar component hosts a dust disk 
with roughly equal amounts of millimeter emission. Neither individual disk is clearly resolved. The 
composite and component-resolved SEDs for this system are shown together in Figure [8j along 
with a potential model for an extincted stellar photosphere from the primary. Since very little 
resolved information is available and the nature of the infrared companion is uncertain, we have 
not modeled the circumstellar material in this system. There is potentially great interest in doing 
so with some modifications to the standard technique presented here. It would be worthwhile to 
better understand the lack of accretion onto the primary despite the presence of a substantial dust 
mass, as well as the potential origin of the very red SED of the companion object. 
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Discussion 



5.1. Constraints on Viscous Evolution 



We have used a radiative transfer modeling technique to extract the dust density structures for 
a significant sample of protoplanetary disks based on high angular resolution observations of their 
millimeter continuum emission. Those data and the SEDs for each source were reproduced well 
using a simple model for the two-dimensional density structure of a viscous accretion disk, with 
a parametric surface density profile that varies with radius like a power-l aw with an exponential 
taper at large radii (see Eq. [1]). This specific form of S corresponds to the iLynden-Bell &: Pringk 
(j 19741 ) similarity solutions for the evolution of a thin accretion disk in Keplerian r otation around a 



(stell ar) point mass, with a viscosity v oc W that does not vary with time (see also lHartmann et al 



19981 ). This type of model - without a sharply truncated outer edge - is favored observationally, 
as it provides a natural explan ation for the optical absorption profiles of silhouette disks in Orion 
(JMcCaughrean fe 0'Defllll996l ) and reconciles the apparent discre pancies in the obser ved spatial 
extents of the dust and CO line emission for nearby resolved disks (JHughes et al.l 1200a ). 



In the context of these models, Figure [3b demonstrates that there is a relatively narrow dis- 
tribution of the parameter 7 that describes the spatial distribution of mass in the disk, consistent 
with a median (7} = 0.9 ± 0.2 when the sample is considered together (estimates for individual 
disks range from 7 = 0.4-1.1). These values fall at the high end of the wider distribution of E 
gradients inferred bv llsella et al.l (|2009l . where 7 ranges from -0.8 to 0.8) using a 1.3 mm continuum 
survey of Taurus disks with slightly poorer angular resolution (0.7-1.0"), and therefore probing the 
mass at larger disk radii. Although there is substantial overlap in the 7 distributions from both 
surveys (particularly when the disks with large central cavities are ignored), the differences can be 
attributed to the distinct approaches for interpreting the data and there is no obvious means of 
reconciliation at present. Focusing on our results, we infer that the mass in ~1 Myr-old disks has a 
similar radial distribution regardless of the wide range of total masses {M<i ~ 0.001-0.1 M ) probed 
in this sample. This implies that whatever mechanism is responsible for generating the viscosity in 
these disks insures that it has a roughly linear dependence on radius (7 ~ 1). 

Since the physical origin of that viscosity is unclear, we adopt a simple prescription v = ac s H, 
where c s is the sound speed and the viscosity co efficient a describes the efficiency of the angular 
momentum transport (IShakura fc Sunyaevlll973l ). Using the density structures and measurements 
of mass accretion rates M* (see Table [6] for references), we can estimate the value of the viscosity 
coefficient a(R) « M*(i?/i? c ) 7 /37rS c c s i7 (see the Appendix in Paper I for details). The resulting 
distribution of a at R = 10 AU is shown in Figure EJ with values for individual disks listed in 
Table Note that for the typical midplane temperature distribution T oc R~ q , the viscosity 
coefficient is proportional to R 1 /c s H ~ R z , where z ~ 7 + q/2 — ip — 1. The model parameters 
and temperature distributions inferred for this sample have z ~ 0.0 ± 0.3, meaning a does not 
vary by more than a factor of ~2-4 across the disk. The inferred distribution of these viscosity 
coefficients appears bimodal, with peaks just below a ~ 0.001 and above 0.01. However, given 
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the small sample, those peaks are not statistically significant; they simply represent the clustering 
of the accretion rates in this sample. The range of inferred a values is in reasonable quantitative 
agreement with magnetohydrodynamics simulations where the vis cosity is generated by turbulence 



from the magnetorotational i nstability in slightly ionized disks (jHawlev et al.l Il995l ; IStone et al 
19961 ; Fleming &: Stond 120031 ; iFromang et al.l 120071 ). However, we should caution that these a 



values are not yet well constrained: they suffer the combined uncertainties in the accretion rates 
and disk structure parameters, as well as from systematic issues with the inherent assumptions in 
their derivation (e.g., the dust traces the gas, viscous heating is insignificant, etc.). They represent 
only a first exploration of an improved empirical understanding of disk viscosities. 

The interactions of gravitational and viscous torques control the evolution of disk structure 
over the vast majority of the disk lifetime. This viscous evolution process has two important, 
observable effects on that structure. First, the coupling of the viscosity and the Keplerian orbital 
shear drives a net mass flow toward small radii, where material can be magnetically channeled 
onto the central star. That mass flow produces bright H emission lines, and the shock generated 
when it impacts the ste llar surface gives rise t o a strong ultraviolet continuum; both tracers ca n 



be used to estimate M* (JMuzerolle et al.lll998al lbl: iGullbring k, CalvetJ[l99^ : lGullbring et aljlioooh . 



To compensate for the angular momentum dissipated in that process, the material in the outer 
disk (beyond some radius Rt] see Table [6] and Paper I) is spread to larger radii. The results of 
that viscous diffusion - decreased average densities (i.e., Mj) and increased sizes (R c ) - can be 
probed with resolved observations of the millimeter continuum emission. In principle, tracking 
these observationa l signatures as a fu nction of time could provide strong constraints on the viscous 
evolution process. Ilsella et al.1 (120091 ) claim a significant correlation between their parameterization 
of R c and stellar age in a sample of Taurus disks. No such correlation is evident in our larger sample 
of Ophiuchus disks. More importantly, the search for such a trend within a single star- forming region 
is probably premature; the ages of individual yo ung stars can not be determined with sufficient 
accuracy to infer such evolutionary behavior (e.g., lHillenbrandll2009l . and references therein), even 
though it may exist. Rather, the inferred disk structures could be considered representative of the 
diversity of evolutionary states and/or initial conditions at a "snapshot" in the evolution sequence 
corresponding to some median cluster age (in this case, ~lMyr). 

While there may not be much difference in the shapes of the £ profiles (7), this sample includes 
disks with a wide variety of masses (Mj) and characteristic size scales (R c )- As was demonstrated in 
Figure [7J the disks with brighter millimeter emission have systematically larger masses and sizes. 
This relationship is illustrated more directly in Figure [TU1 with a significant correlation (3.3 a; 
Spearman rank coefficient of 0.85) between R c and Mj that can be approximated as a power-law, 
Md ex i?J -6±0,3 . The physical origin of the correlation is not clear. One potential explanation is 
that it reflects the range of initial conditions and viscous properties inherited at the disk formation 
epoch. In that case, note that the sense of the correlation is almost perpendicular to the evolutionary 

— 1/2 

paths for individual disks that conserve angular momentum, Md oc R c . The spread along the 
correlation is then representative of the distribution of angular momenta imparted by the collapse of 
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the parent molecular cloud cores. Following llsella et al.l (|2009l ). we can estimate the specific angular 
momenta in these cores that are needed to reproduce the inferred disk structures. Assuming the 
centrifugal radius in the disk corresponds to the radius that contained ~90% of the disk mass at 
the formation epoch (here defined as approximately twice an initial scaling radius, ~2i?i), then the 
specific angular momentum can be written 

0.5 / ,,.r \ -1.5 / n \ 2 



10 



20 



R 



25 AU 



M 



Rr 



0.1 pc 



cm s 



(2) 



where R CO re is the core radius (see iHueso fc Guillotl 120051 ) . Precise values for the initial scaling 
radius are uncertain (see below), but we can estimate upper limits on the j values since R\ < 
R c by definition (see Paper I). Substituting R\ = R c , R core = 0.1 pc, and the stellar masses 
in Table into Eq. (2), we find logj ~ 19.7-20.9 in units of cm 2 s _1 . Radio observations of 
rotating molecular cloud cores have been used to i nfer similar or hig her values, logj ~ 19.6-22.2 



(Goodman et al. 



1993 



Barranco fc Goodmanlll998l ; ICaselli et al.ll2002l ). However, recent numerical 



simulations have de monstrated tha t those measurements tend to overestimate j by roughly an 
order of magnitude (JDib et al.ll2010l ). If that is the case, then the range of disk angular momenta 
inferred here are actually in good quantitative agreement with those determined f or analogs o f thei r 
parental molecular cloud cores. The low end of the j distribution postulated by iDib et al.l (|2010i ) 
and not probed here could then be explained by a selection effect (perhaps our flux-limited sample 
only recovers disks with large j values), an unspecified mechanism for angular momentum loss (e.g., 
outflows or magnetic braking), and/or that R\ <C R c in at least some cases. 

Unfortunately, an unambiguous explanation of the shape of the M$-R c correlation in this 
scenario is not possible, due to the fundamental degeneracy between the initial disk structures 
and the rates at which they evolve. The values of M d and R c change with time such that 
M d = M di0 /T 1/2(2_7) and R c = i?iT 1/(2 ~ 7) , where M dfl and R ± were the initial disk mass and 
characteristic radius, respectively, and T is a dimensionles s parameter that track s how many vis- 
cous timescales (t s ) have elapsed, T = l + t/t s (see Paper h lHartmann et al.lll9 98). That evolution 
is described by 2 equations with 3 unknown parameters, {M^.o, Ri, t s }: the evolution rate can 
not be uniquely disentangled from the initial conditions. However, the lack of significant outliers 
to the correlation may be related to the viscous evolution process itself. The "small and massive" 
region in the upper left of Figure [10] is depopulated because these disks spend only a short time 
with such compact density configurations early in their evolution (t <C lMyr). Likewise, the "large 
and low-mass" region to the lower right in this parameter-space is empty because these ~1 Myr-old 
disks simply have not had enough time to evolve into it. 

Alternatively, the inferred M d -R c correlation could be the result of coupling the viscous evolu- 
tion scenario described above with an internal dissipation process. Recent theoretical calculations 
have suggested that the high-energy (particularly far-ultraviolet, or FUV) irradiation of the disk 
surface by the central star can drive significant mass loss from lar ge disk radii {R > 30 AU) in a 
photoevaporative flow (JGorti fc Hollenbachll2009l ; iGorti et al.ll2009l ). In this scenario, the smaller, 
low-mass disks in this sample should be associated with an advanced stage in this photoevaporation 
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process, perhaps because their central stars have more intense FUV radiation fields. Unfortunately, 
direct constraints on the FUV emission from these extincted sources in the Ophiuchus clouds are 
rare. Instead, we can approximate the FUV luminosity as the sum of a fracti on of the accretio n 



luminosity (oc M*M*/R*) and a chromospheric component («5x 10 -4 L*; see lGorti et al.l 120091 ). 
However, we find no evidence for an expected anti-correlation between disk masses/sizes and the 
FUV (or X-ray) luminosities in this sample: if anything, the larger and more massive disks tend 
to have more intense FUV radiation fields, whereas the less massive, smaller disks show a wide 
range of FUV luminosities. While this certainly does not exclude a role for photoevaporation in 
shaping disk structures, it does suggest that diverse initial conditions and viscous timescales are 
more significant contributors to the observed Md-R c correlation at an age of ~1 Myr. 

Comparisons of disk properties in different stellar populations will be required to definitively 
assess the relative impacts that initial conditions, viscous evolution, and dissipation effects like 
photoevaporation have on shaping disk density distributions over time. The key to disentangling 
the initial conditions from the viscous evolution timescale is to search for and characterize structural 
correlations, like the Md-R c relationship found here, in disks around stars with a range of ages. In 
a scenario where viscous evolution dominates, the disks embedded in the envelopes of Class 0/1 
sources should be massive and compact, providing direct constraints on the distributions of initial 
masses (M^ ) and sizes {R\). Conversely, the disks in older star-forming regions are expected to 
be larger and less massive, unless or until FUV photoevaporation can effectively dissipate the mass 
reservoir in the outer disk. Tracking the evolution of disk masses and sizes in these older stellar 
populations can provide crucial insights into the distribution of viscous timescales, as well as into 
the timeframe where photoevaporation may dissipate much of the disk mass reservoir. To probe 
disk structures in these older star clusters that are typically at larger distances, instruments with 
increased sensitivity to millimeter continuum emission and access to higher angular resolution will 
be required. Fortunately, there is great promise for a breakthrough in our understanding of viscous 
evolution and disk dissipation in the near future, based on extensive millimeter continuum surveys 
of disk populations with the Atacama Large Millimeter Array (ALMA). 



5.2. Implications for Planet Formation 

Ultimately, the drive to understand this viscous evolution process lies in the desire to constrain 
how mass is re-distributed over time, and therefore where and when the conditions in a given disk 
are suitable for planet formation. All planet formation models require that the disk exceeds some 
density threshold as a necessary (but perhaps not sufficient) condition for making a planet. In light 
of the new constraints on disk densities presented here and in Paper I, it is nat ural to com pare 



them with the two basic recipes for gia nt planet formation; disk instability (e.g., lBosslll997l ) and 



core accretion (e.g.. iPollack et al.lll996l ). In the disk instability model, an over-dense region with 



a sufficiently short cooling time can fragment out of the global disk struct ure and precipitate a 



bound gaseous protoplanet very rapidly (~10 3 4 yr; see the recent review bv lDurisen et al.l 120071 ). 
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Typically those conditions are only met at large disk radii (jBolev et al.ll2006l ); such protoplanets 
must then migrate inward to reproduce the shorter period orbits observed in the Solar System 
and elsewhere. In the core accretion model, a giant planet is produced from a relatively slow (~1- 
10 Myr) collisional growth proces s that first builds up a large solid core a nd then rapidly accretes 
a massive gaseous envelope (e.g.. iHubickvi et al.ll2005l : lAlibert et al.ll2005l ). 



Assuming the dust traces 1% of the gas mass, the disk structures presented here are stable 
against gravitational fragmentation, with the minimum local Toomre Q yalues ranging from ~5-50 
in most cases (Q = CsQ/ttGE, where f2 is the Keplerian orbital velocity; lToomrdll964l ). The large, 
massive disks around DoAr 25 and GSS 39 are potentially exceptions, with Q values approaching 
2 at radii of ~65 and 150 AU, respectively. However, we should caution that the disk structures 
in both cases are particularly uncertain (see §4.2 of Paper I). While none of the sample disks are 
clear candidates for the efficient operation of the disk instability mechanism for planet formation in 
the current epoch (i.e., Q < 1.7), their structures may have been le ss stable at earlier times. Even 
without invoking some kind of mass loading mechanism (|Boleyll2009l ) , viscous evolution implies that 
these disks originally had denser, more compact structures, which might suggest a formerly lower 
minimum Q-value and/or radius where the disk is least stable to fragmentation. A quantitative 
exploration of the history of Q(R) in these disks grounded in our constraints on their viscous 
properties may well be worthwhile, but would require new radiative transfer calculations that 
incorporate some model for the thermal evolution of the disk material in each case. 

A similar comparison with the core accretion model is more challenging, primarily because 
resolution limitations do not yet permit a direct observation of the disk material inside R ~ 20 AU. 
Typically, core accretion model calculations impose a scaled Minimum Mass Solar Nebula (MMSN) 
density structure for an initial disk, where S oc R~ 3 < 2 ( Weidenschillinall977l ) and a total mass set to 
produce a formation efficiency commensurate with the observational constraints on disk lifetimes. 
Although the MMSN surface density profile is too steep compared to those derived here (where 
S oc i?™ 1 at the relevant radii), Figure [6] demonstrates that the surface densities for the sample are 
generally compatible with those expected in the 20-40 AU region of the primordial Solar disk. The 
shape of the S profile should not adversely affect the likelihood of planet formation if sufficient mass 
is available, but it is expected to ha ve an impact on the orbital architecture and mig r ation properties 



of any resulting p lanetary system (IKokubo fc Idall2002l : IChambers k, Cassenll2002l : iRaymond et al 



20051 : ICridall2009l ). Some core accretion sim ulations have started to employ viscous d isk density 
structures similar to those derived here (e.g., lAlibert et alj|2005l : iHueso fc Guillotll2005l ). but they 
still must populate that density structure with solid bodies orders of magnitude larger than those 
responsible for the observed millimeter emission. 

This last point highlights an important uncertainty in deriving densities from observations of 
young circumstellar disks: the conversion of emission to mass does not account for particles much 
larger than the observing wavelength. Since this potential for substantially under-estimating the 
disk densities was already elaborated in Paper I, we will not dwell on it here. The problem is not 
likely to be severe in the regions that can currently be probed observationally (R > 20 AU), as 
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even optimistic estimates of pa rticle growth timescales at th ose radii do not predict a substantial 
population of large solids (e.g., iDullemond &: Dominikll2005l ), However, at smaller disk radii those 
growth timescales can be considerably shorter, generating a dust emissivity that varies with radius 
inside the planet formation zone. More generally, it is worthwhile to emphasize that the observations 
are not yet sensitive to the midplane material for R < 20 AU. Inside that resolution limit, the 
densities could be dramatically different than inferred here if, for example, the dust size distribution 



is changed or the effective viscosity profile is modified (e.g. JZhu et al.ll2010l ). New data from ALMA 
and the Expanded Very Large Array (EVLA) will help alleviate these uncertainties at small radii, 
and should also help empirically guide our understanding of disk viscosities. 

But despite our general ignorance of the inner disk structures, a subset of young disks show 
unambiguous evidence for evolution in the planet formation zone. In this SMA survey of 17 disks, 
selected primarily for their millimeter luminosities, 4 of them - around SR 24 S, SR 21, WSB 
60, and DoAr 44 - have large central cavities with significantly diminished dust emission. The 
characteristic ring-like emission morphologies and distinctive visibility nulls for these disks indicate 
that the radial transition in the millimeter optical depth across the cavity edge is rather sharp, 
commensurate with our simple models that employ a density contrast of ~10 2 -10 4 within a small 
radial range. In 3 of these 4 disks, a small amount of micron-sized dust grains must reside inside the 
cavity to account for their observed infrared excess emission. Outside of their central cavities, these 
disks have £ profiles comparable to the disks with continuous dust distributions. At this point, 
the physical mechanism re sponsible for the dust evolution in these disk cavities remains a subject 
of active debate (e.g., see iD'Alessio et alJl200a ; iNajita et al.1 120071 ). They could be produced by 
an abrupt emissivity decre ase due to grain growth in a localized region of low turbulence (e.g., an 
MRI-inactive "dead" zone: ICieslal 120071 : IZhu et al.ll2010l ). Alternatively, they may represent a true 
lack of material in the inner disk, cleared out by a photoevaporative wind and viscous draining 



( Alexander Sz Armitagd 200711 or tidal interactio ns with faint companions - perhap 



planetary systems (e.g.. lLubow &: D , Angeldl2006l ) 



s even young 



6. Summary 

We have conducted a high angular resolution (down to FWHM scales of / .'3 ~ 40 AU) Sub- 
millimeter Array survey of the 880 fim continuum emission from 17 p rotoplanetary disks in the 
~1 Myr-old Ophiuchus star- forming region (see also Andrews et al.ll2009l ). Using a two-dimensional 
parametric model for the structure of a viscous accretion disk and a Monte Carlo radiative transfer 
code to simultaneously reproduce the millimeter visibilities and broadband SEDs, we have measured 
the spatial distribution of mass in these disks. The key conclusions from this survey include: 



1. Assuming the viscosity scales with radius like a power-law (y oc W) and the surface densities 



like £ oc (R/R c ) 7 exp [— {R/R c ) 7 ], the disks in this survey exhibit a relatively narrow 
range of density /viscosity gradients, 7 = 0.4-1.1. Taken together as a sample, these results 
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are consistent with a single (median) value, (7) = 0.9 ± 0.2, independent of the millimeter 
luminosity (disk mass) or the stellar properties. 

2. The disk masses (M^) and characteristic radii (R c ) are correlated: massive disks tend to be 
larger (i.e., brighter disks are more resolved). The 3.3 a correlation is described well by a 
power-law relation, Mj oc i?J ,6±0 ' 3 (with a residual scatter of ~0.3dex). Rough estimates 
of the range of angular momenta in these disks are found to be comparable to the values 
inferred for rotating molecular cloud cores. If it can be measured in disk samples with a 
well-established range of ages, the Md-R c relationship can be used to deduce a mean viscous 
timescale and the range of initial conditions imparted by the disk formation process. 

3. Based on the derived physical conditions in these disks and literature measurements of 
their mass accretion rates , we have made crude estimates of their viscosity coefficients (a; 



Shakura &: Sunvaevl 119 73). The inferred a values range from 0.0005-0.08, with a median 
(a) ~ 0.01. This range of values is commensurate with numerical simulations where the effec- 
tive viscosities are generated by turbulence from the magnetorotational instability in slightly 
ionized disks. 

4. The inferred disk surface densities in the ~20-40 AU range are in good quantitative agreement 
with those expected for the primordial disk that gave rise to the giant planets Uranus and 
Neptune (i.e., the outer part of the Minimum Mass Solar Nebula). Current limitations on 
angular resolution prevent more direct constraints on E in the inner disk (R < 20 AU) . All 
of the sample disks appear to be gravitationally stable (typical minimum Toomre Q values 
range from ~5-50). However, these results do not rule out periods where the disk instability 
mechanism for planet formation at earlier epochs (<ClMyr). 

5. Regions of significantly diminished millimeter emission were resolved at the centers of 4/17 
disks in the sample (SR 24 S, SR 21, DoAr 44, and WSB 60). Simple models for these 
"transition" disks can reproduce the observations well if the densities are decreased by a 
factor of >100 inside a radius R csv ~ 20-40 AU. The disk properties outside these central 
cavities are comparable to the other disks in the sample. The significant infrared excesses 
still present in 3/4 of these disks with resolved central cavities indicate that some small, warm, 
dust particles remain near their central stars. The implied dust structures in those cases are 
commensurate with the opening of large gaps in the planet formation zones of these disks. 
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Table 1. SMA Observing Journal 



Name 


a [J2000] 


5 [J2000] 


Array 


UT Date 


Alt. Name 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Elias 24 


16 26 24.08 


-24 16 13.7 


V 

c 


2009 March 14 
2009 May 2 


WSB 31, YLW 32 


SR24S 


16 26 58.51 


-24 45 37.0 


V 

c 


2009 March 25 
2009 May 4 


DoAr 29, Haro 1-7 


SR4 


16 25 56.16 


-24 20 48.5 


V 

V 

c 


2009 March 12 
2009 March 25 
2009 May 4 


V2058 Oph, AS 206 


SR 13 


16 28 45.27 


-24 28 19.2 


V 

c 


2009 March 12 
2009 May 2 


V853 Oph, HBC 266 


WSB 52 


16 27 39.44 


-24 39 15.7 


V 
V 

c 


2009 March 2 
2009 March 28 
2009 May 4 


ROXs 27 


DoAr 33 


16 27 39.01 


-23 58 18.9 


V 
V 

c 


2009 March 2 
2009 March 28 
2009 May 4 


WSB 53, ROXs 30C 


WL 18 


16 26 48.98 


-24 38 25.4 


V 

c 


2009 March 13 
2009 May 2 


GY 129 


DoAr 24 E 


16 26 23.37 


-24 20 59.8 


V 

c 


2009 March 13 
2009 May 2 


Elias 22, GSS 31 



Note. - Col. (1): Disk name. Cols. (2) & (3): Centroid coordinates, determined as 
described in the text (§2). The coordinates listed for DoAr 24 E correspond to the optically- 
visible component to the northwest. Col. (4): Array configuration; V = very extended (68- 
509 m baselines) and C = compact (16-70 m baselines). Col. (5): UT date of observation. 
Col. (6): Common alternative identifications. 
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Table 2. Continuum and CO Synthesized Map Properties 



Disk 




continuum 






CO J=3-2 






F v 


e b 


PA b 


rms 


o b 


PA b 




[mJy] 


["} 


[°] 


[Jy] 


["] 


[°] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Elias 24 


890 ±3 


0.65 x 0.51 


18 


0.19 


2.13 x 1.42 


51 


SR24S 


545 ±3 


0.51 x 0.43 a 


89 


0.25 


2.15 x 1.44 


52 


SR4 


150 ±2 


0.50 x 0.43 


15 


0.24 


2.12 x 1.45 


53 


SR13 


149 ±3 


0.63 x 0.52 


39 


0.19 


2.15 x 1.41 


50 


WSB 52 


147 ±3 


0.50 x 0.43 


52 


0.25 


2.16 x 1.44 


52 


DoAr 33 


80 ±2 


0.52 x 0.43 


67 


0.25 


2.17 x 1.42 


53 


WL 18 


51 ±3 


0.61 x 0.52 


55 


0.19 


2.19 x 1.39 


50 


DoAr 24 E 


49±2 b 


0.60 x 0.54 


58 


0.19 


2.21 x 1.37 


50 



Note. ■ - Col. (1): Disk name. Col. (2): Integrated continuum flux 
density and rms noise level per beam (does not include the ~10% flux 
calibration uncertainty). Cols. (3) &: (4): FWHM dimensions and po- 
sition angle (measured east of north) of the synthesized beam for the 
continuum maps shown in Figure [2j Col. (5): The rms noise level per 
beam for an individual 0.70 km s _1 channel in the synthesized CO J=3— 2 
channel maps. Cols. (6) &: (7): FWHM dimensions and position angle of 
the synthesized beam for the channel maps. 

a The high-resolution inset map of the SR 24 S disk shown in Figure 
[2] was generated with uniform visibility weighting and has a synthesized 
beam with dimensions 0'.'37 x 0726 at PA = 14°. 

The infrared companion source to the southeast is just slightly fainter, 
with an integrated flux density of 34 ± 2 mJy. 
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Table 3. Stellar Properties 



Name 


SpT 


Ay 


Teg 


R* 


Lit 


Af, 


age 


ref 






[mag] 


[K] 


[Rq] 


[L Q ] 


[M ] 


[Myr] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


AS 205 


K5 


2.9 


4250 


3.7 


4.0 


1.0 


0.5 


1 


Elias 24 


K5 


8.7 


4250 


4.2 


5.1 


1.0 


0.4 


2 


GSS 39 


MO 


15 


3850 


2.3 


1.0 


0.6 


1.0 


3 


AS 209 


K5 


0.9 


4250 


2.3 


1.5 


0.9 


1.6 


4 


DoAr 25 


K5 


2.7 


4250 


1.7 


0.8 


1.0 


3.8 


2 


SR24S 


K2 


7.0 


4990 


2.8 


4.4 


2.0 


2.4 


3 


WaOph 6 


K6 


3.6 


4205 


3.2 


2.9 


0.9 


0.7 


5 


SR 21 


G3 


6.3 


5800 


3.3 


11 


2.0 


4.7 


1 


VSSG 1 


MO 


14 


3850 


3.1 


1.9 


0.6 


0.7 


6 


WSB 60 


M4 


3.5 


3370 


1.3 


0.2 


0.3 


3.0 


2 


DoAr 44 


K3 


2.3 


4730 


1.7 


1.3 


1.4 


7.1 


7 


SR4 


K7 


1.3 


4060 


2.0 


1.0 


0.7 


1.8 


3 


SR13 


M4 


0.0 


3370 


1.9 


0.4 


0.3 


1.5 


2 


WSB 52 


Ml 


5.0 


3750 


1.9 


0.6 


0.5 


1.5 


3 


DoAr 33 


K4 


3.7 


4470 


1.8 


1.2 


1.3 


4.4 


7 


WL 18 


K7 


11 


4060 


1.1 


0.3 


0.8 


11 


2 



Note. - Col. (1): Disk name (those in italics were modeled in 
Paper I). Col. (2): Spectral type. Col. (3): Extinction. Col. (4): 
Effective temperature. Col. (5): Radius. Col. (6): Luminosity . 
Col. (7) and (8): Mass and age estimated from the lSiess et al.l ()2000l ) 
pre-main-seque nce models. Col. (9): Li terature refere n ces fo r SpT 



and Ay. \1] - IPrato et al.l (I2003ll. [21 - IWilking et all (I2005| N I. [31 



Luhman k Ri ekc (1999|), [41 -iHerbig k Belli (11988ft. [51 -Eisner et al 



(J2005I ). [6] - iNattaet al.1 (J2006I ). [7] - Bouvier k Appenzellerl (J1992 ) 
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Table 4. Disk Structure Model Parameters: Continuous Cases 



Name 


M d 


7 


Re 


#100 


v> 


R\n 


i 


PA 


x 2 - 

AyiS 


A.scd 




[M ] 




[AU] 


[AU] 




[AU] 


[°] 


[°] 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


AS 205 


0.029 


0.9 


46 


19.6 


0.11 


0.14 


25 


165 


2.1 


3.7 


Elias 24 


0.117 


0.9 


127 


8.6 


0.03 


0.16 


24 


50 


2.0 


1.5 


GSS 39 


0.143 


0.7 


198 


7.3 


0.08 


0.07 


60 


110 


1.9 


32 


AS 209 


0.028 


0.4 


126 


13.3 


0.10 


0.09 


38 


86 


1.7 


2.4 


DoAr 25 


0.136 


0.9 


80 


6.7 


0.15 


0.06 


59 


112 


1.9 


9.2 


WaOph 6 


0.077 


1.0 


153 


4.4 


0.06 


0.12 


39 


171 


1.8 


1.8 


VSSG 1 


0.029 


0.8 


33 


9.7 


0.08 


0.10 


53 


165 


1.8 


12 


SR4 


0.004 


0.8 


20 


19.9 


0.23 


0.07 


50 


39 


1.9 


8.8 


SR 13 


0.012 


1.0 


26 


12.6 


0.07 


0.04 


32 


42 


2.0 


1.9 


WSB 52 


0.007 


1.1 


26 


14.5 


0.19 


0.06 


46 


120 


1.9 


6.6 


DoAr 33 


0.007 


1.1 


38 


2.7 


0.06 


0.07 


43 


102 


1.9 


22 


WL 18 


0.011 


0.8 


14 


5.6 


0.07 


0.04 


48 


115 


1.9 


13 



Note. - - Col. (1): Disk name (those in italics were modeled in Paper I). 
Col. (2): Disk mass assuming a 100:1 gas-to-dust mass ratio. Col. (3): Ra- 
dial surface density gradient. Col. (4): Characteristic radius. Col. (5): Scale 
height at 100 AU. Col. (6): Radial scale height gradient. Col. (7): Fixed inner 
radius. Col. (8): Fixed inclination. Col. (9): Fixed major axis position angle. 
Col. (10): Reduced % 2 statistic comparing the model fit with the continuum 
visibilities alone. Col. (11): Same as Col. (10), but for the SED alone. 
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Table 5. Disk Structure Model Parameters: Central Cavity Cases 



Name 


M d 


7 


R r 


#100 


$ 


-^cav 


^cav 


i 


PA 


X 2 - 

A. vis 


-^sed 




[M e ] 




[AU] 


[AU] 




[AU] 




[°] 


[°] 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


SR24S 


0.042 


0.8 


40 


4.0 


0.02 


32 


0.0001 


57 


31 


2.1 


11 


SR 21 


0.005 


0.9 


17 


7.7 


0.26 


37 


0.005 


22 


110 


1.7 


7.2 


WSB 60 


0.021 


0.8 


31 


11.0 


0.13 


20 


0.01 


25 


117 


1.8 


3.0 


DoAr 44 


0.017 


1.0 


80 


3.5 


0.04 


33 


0.0001 


45 


75 


1.8 





Note. -- Cols. (1-6): Same as for Table HI Col. (7): The cavity radius, marking the 
outer edge of the diminished inner disk densities (see Paper I or the discussion on SR 
24 S in §4.2). Col. (8): The density reduction scaling factor inside the radius i? ca v 
Cols. (9-12): Same as for Tabled cols. (8-11). 
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Table 6. Viscous Disk Properties 



Name 


[M yr" 1 ] 


Rt 

[AU] 


a 


ref 


(1) 


(2) 


(3) 


(4) 


(5) 


AS 205 


8 x 1(T 8 


23 


0.005 


1 


Elias 24 


2 x 1(T 7 


62 


0.03 


2 


GSS 39 


7 x 1(T 8 


95 


0.03 


2 


AS 209 


9 x 1(T 8 


61 


0.08 


3 


DoAr 25 


3 x 1(T 9 


39 


0.0005 


4 


WaOph 6 


1 x 1(T 7 


78 


0.05 


5 


VSSG 1 


1 x 1(T 7 


16 


0.02 


2 


SR4 


6 x 1(T 8 


10 


0.02 


2 


SR13 


3 x 1(T 9 


13 


0.0005 


2 


WSB 52 


4 x 1(T 9 


13 


0.001 


2 


DoAr 33 


3 x 1(T 10 


20 


0.0006 


6 


WL 18 


9 x 1(T 9 


7 


0.003 


2 


SR24S 


3 x 10" 8 


19 


(0.006) 


2 


SR 21 


< 2 x 1(T 9 


9 


(<0.0009) 


2 


WSB 60 


1 x 1(T 9 


15 


(0.0002) 


2 


DoAr 44 


9 x 1(T 9 


40 


(0.01) 


7 



Note. - Col. (1): Disk name (those in italics 
were modeled in Paper I). Col. (2): Accretion rate. 
Col. (3): Radius of mass flow reversal. Col. (4): 
Vis cosity coeffic i ent. Col. ( 5): References for M* : 



,Tohns-Kru 1 et al.1 (l200(t l [4] - JLuhm 




Pratoet alJteOO.Tl. \2) - JNatta et a l. (200 6J), [31 



Eisner et al.1 (I2005l'l [6] 



[7] - lEspaillat et al.1 (|2010|). 



an &: Rieke 



Cieza et al. 
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Fig. 1. — The millimeter continuum flux densities and visual extinctions for young stars in the 
Ophiuchus molecular clouds (Andrews & Williams 2007a, b; Cieza et al. 2010; additional Ay esti- 
mates from Furlan et al. 2009). The targets selected for this sample are shown in red (880 /um flux 
densities are from Tableland Paper I), and the regions that do not meet the selection criteria are 
shaded grey. Blue arrows represent 3<r upper limits and black points mark sources that have not 
been observed with the SMA at high angular resolution. Embedded sources with uncertain spectral 
classifications are represented with a lower limit of Ay = 50. When 880 ^m flux densities are not 
available, 1.3mm data were scaled up by a factor of ~2.2 (see text); the right-hand ordinate axis 
represents an approximate conversion of the flux density scale to total disk masses (see Andrews & 
Williams 2005; 2007b). 
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Fig. 2. — Aperture synthesis images of the 880 ^m continuum emission from the 8 disk targets in 
the expanded sample. Each panel is 4" (500 AU) on a side. Contours start at 3 a and increase in 
3cr intervales (5cr for Elias 24 and SR 24 S only; rms values are given in Table [2]). The synthesized 
beams are shown in the lower left corner of each panel. Note the detection of roughly equal emission 
levels around both the optically- visible primary star DoAr 24 E and its infrared companion to the 
southeast, as well as the asymmetric ring-like emission morphology for the disk around SR 24 S 
(shown separately in detail). 
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Fig. 3. — The derived distributions of the disk structure parameters for the composite sample 
(combining the results presented here and in Paper I). From left to right are the disk masses (M^), 
radial surface density gradients (7), characteristic radii (R c ), scale- heights at 100 AU (Hiqq), and 
the radial scale-height gradients (ip). The contributions of the four disks with diminished millimeter 
emission in their central regions are hatched. 
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Fig. 4. — Comparison of the data with the best-fit disk structure models. The left panels show the 
SMA continuum image, corresponding model image, and imaged residuals (data— model). Contours 
are drawn at the same 3<r intervals in each panel. Crosshairs mark the disk centers and major axis 
position angle; their relative lengths represent the disk inclination. The right panels show the 
broadband SEDs and deprojected visibility profiles, with best-fit models overlaid in red. The input 
stellar photospheres are shown as blue dashed curves. 
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Fig. 5. — Same as Figured] for the SR 24 S disk, which features a resolved central emission cavity. 
The modifications to the standard disk model made for this source in particular are detailed in 
§4.2. The inset images are to scale, and were synthesized by uniformly weighting the visibilities to 
highlight the structure at higher angular resolution. 
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continuous disks 




100 
R [AU] 



1000 



Fig. 6. — Surface density profiles for the sample, where the targets with continuous emission 
distributions (top) and central emission cavities (bottom) are separated for clarity. The new model 
results presented here are highlighted in color, while those from Paper I are shown in black. The 
light grey band out to a radius of 20 AU marks the resolution limit of the survey. Dark gray 
boxes represent the surface de nsities extrapolated fo r Saturn, Uranus, and Neptune in the standard 
Minimum Mass Solar Nebula (|Weidenschillinglll977l ) . The density profiles inferred from dust tracers 
have been scaled up assuming a gas-to-dust mass ratio of 100:1. 
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Fig. 7. — A direct comparison of the sample-averaged visibility profiles for the subsets of disks 
with integrated 880 /im flux densities fainter (blue) and brighter (red) than 0.5 Jy. Each profile was 
constructed by averaging the visibilities for individual disks in annular bins after their deprojection 
according to viewing geometry and normalization by their total flux densities. The disks with large 
central emission cavities were excluded for clarity. Error bars correspond to the standard deviation 
in each bin. The significantly larger amount of correlated emission on all spatial scales for the 
fainter subset of disks demonstrates that they are less well resolved, and therefore substantially 
smaller than their brighter counterparts. This empirical relationship is reinforced by the radiative 
transfer modeling, which implies a modest correlation between Mj and R c (see §5.1). 
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Fig. 8. — The SED for the combined DoAr 24 E system (black circles) is shown together with the 
component-resolved SEDs of the optical primary (red) and infrared companion (blue) that lie to 
the northwest and southeast in Figure [H respectively. A model stellar photosphere for the optical 
primary is shown as a dashed red curve. Unlike those shown in Figure HI these SEDs have not been 
corrected for extinction due to the uncertainty in the nature of the infrared companion source. 
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Fig. 9. — The distribution of viscosity coefficients, a, inferred for this sample from the disk struc- 
tures and mass accretion rates. The individual a values listed in Table [6] are still quite uncertain, 
although the broad range is commensurate with simulations where the disk viscosities are sus- 
tained by MHD turbulence. Given the small sample size, a uncertainties, and potential M* bias, 
the bimodal appearance of the distribution is not significant. 
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Fig. 10. — A 3.3(7 correlation between the disk masses and characteristic radii. The solid line 
shows a power-law that provides a reasonable match to this relationship, where M^ ex R^ with 
x = 1.6 ± 0.3. Diamond symbols mark the disks with large central emission cavities, and the cross 
to the upper left is representative of the typical parameter uncertainties. The arrow denotes the 
generic direction of evolution implied by the conservation of angular momentum. The range along 
the correlation is representative of the range of angular momenta incorporated into the disks when 
they formed, but the shape of the correlation encodes both these initial conditions and potentially 
a range of viscous timescales. 



-35- 

REFERENCES 

Alexander, R. D., & Armitage, P. J. 2007, MNRAS, 375, 500 

Alibert, Y., Mordasini, C, Benz, W., & Winisdoerffer, C. 2005, A&A, 434, 343 

Andre, P., & Montmerle, T. 1994, ApJ, 420, 837 

Andrews, S. M., k Williams, J. P. 2005, ApJ, 619, L175 (2005a) 

- 2005, ApJ, 635, 1134 (2005b) 

- 2007, ApJ, 659, 705 (2007a) 

- 2007, ApJ, 671, 1800 (2007b) 

Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C, & Dullemond, C. P. 2009, ApJ, 700, 1502 
(Paper I) 

Barranco, J. A., & Goodman, A. A. 1998, ApJ, 504, 207 

Barsony, M., Ressler, M. E., & Marsh, K. A. 2005, ApJ, 630, 381 

Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Glisten, R. 1990, AJ, 99, 924 

Boley, A. C, Mejia, A. C, Durisen, R. H., Cai, K., Pickett, M. K., & D'Alessio, P. 2006, ApJ, 651, 
517 

Boley, A. C. 2009, ApJ, 695, L53 

Bontemps, S., et al. 2001, A&A, 372, 173 

Boss, A. P. 1997, Science, 276, 1836 

Bouvier, J., & Appenzeller, I. 1992, A&AS, 92, 481 

Brown, J. M., Blake, G. A., Qi, C, Dullemond, C. P., & Wilner, D. J. 2008, ApJ, 675, L109 

Brown, J. M., Blake, G. A., Qi, C, Dullemond, C. P., Wilner, D. J., & Williams, J. P. 2009, ApJ, 
704, 496 

Caselli, P., Benson, P. J., Myers, P. C, & Tafalla, M. 2002, ApJ, 572, 238 

Chambers, J. E., & Cassen, P. 2002, Meteoritics & Planetary Science, 37, 1523 

Chiang, E. I., & Goldreich, P. 1997, ApJ, 490, 368 

Ciesla, F. J. 2007, ApJ, 654, L159 

Cieza, L. A., et al. 2010, ApJ, 712, 925 



-36- 

Crida, A. 2009, ApJ, 698, 606 

Cutri, R. M., et al. 2003, 2MASS All-Sky Point Source Catalog (Pasadena: IPAC) 

D'Alessio, P. et al. 2005, ApJ, 621, 461 

Dartois, E., Dutrey, A., & Guilloteau, S. 2003, A&A, 399, 773 

de Geus, E. J., de Zeeuw, P. T., & Lub, J. 1989, A&A, 216, 44 

Dib, S., Hennebelle, P., Pineda, J. E., Csengeri, T., Bontemps, S., Audit, E., &; Goodman, A. A. 
2010, ApJ, in press (arXiv:1003.5118) 

Draine, B. T., & Lee, H. M. 1984, ApJ, 285, 89 

Dullemond, C. P., & Dominik, C. 2004, A&A, 417, 159 

- 2005, A&A, 434, 971 

Durisen, R. H., Boss, A. P., Mayer, R. L., Nelson, A. F., Quinn, T., & Rice, W. K. M. 2007, in 
Protostars & Planets V, eds. B. Reipurth, D. Jewitt, &: K. Keil (Tucson: Univ. Arizona 
Press), 607 

Eisner, J. A., Hillenbrand, L. A., White, R. J., Akeson, R. L., & Sargent, A. I. 2005, ApJ, 623, 952 

Espaillat, C., et al. 2007, ApJ, 670, L135 

- 2010, ApJ, in press flarXiv: 1005.2365 p 
Evans, N. J., et al. 2003, PASP, 115, 965 
Fleming, T., & Stone, J. M. 2003, ApJ, 585, 908 

Fromang, S., Papaloizou, J., Lesur, G., & Heinemann, T. 2007, A&A, 476, 1123 

Furlan, E., et al. 2009, ApJ, 703, 1964 

Garaud, P. 2007, ApJ, 671, 2091 

Ghez, A. M., Neugebauer, G., & Matthews, K. 1993, AJ, 106, 2005 

Goodman, A. A., Benson, P., Fuller, G., & Myers, P. 1993, ApJ, 406, 528 

Gorti, U., k Hollenbach, D. 2009, ApJ, 690, 1539 

Gorti, U., Dullemond, C. P., & Hollenbach, D. 2009, ApJ, 705, 1237 

Gullbring, E., & Calvet, N. 1998, ApJ, 509, 802 

Gullbring, E., Calvet, N., Muzerolle, J., & Hartmann, L. 2000, ApJ, 544, 927 



-37- 

Hartmann, L., Calvet, N., Gullbring, E., k D'Alessio, P. 1998, ApJ, 495, 385 

Hawley, J. F., Gammie, C. F., k Balbus, S. A. 1995, ApJ, 440, 742 

Herbig, G. H., k Bell, K. R. 1988, in Third Catalog of Emission Line Stars of the Orion Population 
(Santa Cruz: Lick Obs.) 

Herbst, W., Herbst, D. K., k Grossman, E. J. 1994, AJ, 108, 1906 

Hillenbrand, L. A. 2009, in IAU Symp. 258, The Age of Stars, ed. E. Mamajek, et al. (Cambridge: 
Cambridge University Press), 81 

Ho, P. T. P., Moran, J. M., k Lo, K. Y. 2004, ApJ, 616, LI 

Hubickyj, O., Bodenheimer, P., k Lissauer, J. J. 2005, Icarus, 179, 415 

Hueso, R., k Guillot, T. 2005, A&A, 442, 703 

Hughes, A. M., Wilner, D. J., Calvet, N., D'Alessio, P., Claussen, M. J., k Hogerheijde, M. R. 
2007, ApJ, 664, 536 

Hughes, A. M., Wilner, D. J., Qi, C, k Hogerheijde, M. R. 2008, ApJ, 678, 1119 

Hughes, A. M., et al. 2009, ApJ, 698, 131 

Ida, S., k Lin, D. N. C. 2004, ApJ, 604, 388 

Isella, A., Carpenter, J. M., k Sargent, A. I. 2009, ApJ, 701, 260 

- 2010, ApJ, 714, 1746 

Jensen, E. L. N., k Mathieu, R. D. 1997, AJ, 114, 301 

Johns-Krull, C. M., Valenti, J. A., k Linsky, J. L. 2000, ApJ, 539, 815 

Kitamura, Y., Momose, M., Yokogawa, S., Kawabe, R., Tamura, M., k Ida, S. 2002, ApJ, 581, 357 

Kokubo, E., k Ida, S. 2002, ApJ, 581, 666 

Knude, J., k H0g, E. 1998, A&A, 338, 897 

Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodriguez, L. F. 2008, ApJ, 675, L29 

Lombardi, M., Lada, C. J., k Alves, J. 2008, A&A, 480, L785 

Lubow, S. H., k D'Angelo, G. 2006, ApJ, 641, 526 

Luhman, K. L., k Rieke, G. H. 1999, ApJ, 525, 440 

Lynden-Bell, D., k Pringle, J. E. 1974, MNRAS, 168, 603 



38 



Mayama, S., et al. 2010, Science, 327, 306 

McCabe, C, Ghez, A. M, Prato, L., Duchene, G., Fischer, R. S., & Telesco, C. 2006, ApJ, 636, 
932 

McCaughrean, M. J., & O'Dell, C. R. 1996, AJ, 111, 1977 

Mordasini, C., Alibert, Y., k Benz, W. 2009, A&A, 501, 1139 

Motte, F., Andre, P., & Neri, R. 1998, A&A, 336, 150 

Muzerolle, J., Calvet, N., & Hartmann, L. 1998, ApJ, 492, 743 

Muzerolle, J., Hartmann, L., & Calvet, N. 1998, AJ, 116, 2965 (1998b) 

Natta, A., Testi, L., & Randich, S. 2006, A&A, 452, 245 

Najita, J. R., Strom, S. E., & Muzerolle, J. 2007, MNRAS, 378, 369 

Niirnberger, D., Brandner, W., Yorke, H. W., & Zinnecker, H. 1998, A&A, 330, 549 

Padgett, D. C., et al. 2008, ApJ, 672, 1013 

Patience, J., Akeson, R. L., & Jensen, E. L. N. 2008, ApJ, 677, 616 

Pietu, V., Dutrey, A., Guilloteau, S., Chapillon, E., & Pety, J. 2006, A&A, 460, L43 

Pietu, V., Dutrey, A., & Guilloteau, S. 2007, A&A, 467, 163 

Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., Podolak, M., & Greenzweig, Y. 1996, 
Icarus, 124, 62 

Prato, L., Greene, T. P., & Simon, M. 2003, ApJ, 584, 853 

Raymond, S. N., Quinn, T., & Lunine, J. I. 2005, ApJ, 632, 670 

Reipurth, B., & Zinnecker, H. 1993, A&A, 278, 81 

Schaefer, G. H., Simon, M., Beck, T. L., Nelan, E., & Prato, L. 2006, AJ, 132, 2618 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 

Simon, M., et al. 1995, ApJ, 443, 625 

Stanke, T., Smith, M. D., Gredel, R., & Kanzadyan, T. 2006, A&A, 447, 609 

Stone, J. M., Hawley, J. F., Gammie, C. F., & Balbus, S. A. 1996, ApJ, 463, 656 



-39- 

Toomre, A. 1964, ApJ, 139, 1217 

Weidenschilling, S. J. 1977, Ap&SS, 51, 153 

Weingartner, J. C, & Draine, B. T. 2001, ApJ, 548, 296 

Wilking, B. A., Meyer, M. R., Robinson, J. G., & Greene, T. P. 2005, AJ, 130, 1733 

Zhu, Z., Hartmann, L., & Gammie, C. 2010, ApJ, 713, 1143 



This preprint was prepared with the AAS IATpjX macros v5.2. 



